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Here we consider the Husimi function P for the squeezed states and calculate the marginal and 
correlation distribution functions when P is projected onto the photon number states. According to 
/-»/ ' the value of the squeezing parameter one verifies the occurence of oscillations and beats as already 

Qv. , appointed in the literature. We verify that these phenomena are entirely contained in the correlation 

O^ ' function. In particular, we show that since the Husimi and its marginal distribution functions satisfy 

partial differential equations where the squeeze parameter plays the role of time, the solutions (the 
squeezed functions obtained from "initial" unsqueezed functions) can be expressed by means of 
O ' kernels responsible for the "propagation" of squeezing. From the calculational point of view, this 

^^ I method presents advantages for calculating the marginal distribution functions (compared to a direct 

integration over one of the two phase-space variables of P) since one can use the symmetry properties 
of the differential equations. 
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O . I. INTRODUCTION 

o\ 

Q>^ ■ In the end of the eighties, the oscillations of the photon distribution function (DF) of high energy squeezed and 

correlated states were discovered in MM ; the authors of |l| studied the oscillatory behavior of the single-mode squeezed- 
state Husimi function (HP) projected into photon number states, P„(p, g; A, 0) = \{n\p,q; X,(f>)\'^, where p and q are 
the space variables associated with the two quadratures of a monochromatic electromagnetic (EM) field, A is the 
r^ ■ squeezing parameter and is a rotation angle in phase space. They suggested that ioi p = 0, q — 7\/2, cj) — and a 
C^ ' fixed A = 21, such oscillatory behavior can be explained in terms of quantum interference effects and were taken as a 
^ : signature of a nonclassical state. More recently, Gagen [gj generalized this study to incorporate interference structures 
^ ' in the Bohr-Sommerfeld trajectories associated with a superposition of quantum states. On the other hand, Dutta et 
^ , al. [H verified an additional feature present in Pn{p, q', A, 0), when plotted as function of n: this distribution exhibits 
JJ"! ' the structure of beats (collapses and revivals) at large values of n (> 10) for A = 201, p^ + q^ = 98 and w 7r/2. 
rN Compared to the value A used in IT], the high value for A is crucial for the occurence of beats. These oscillations were 

jrt attributed to interference effects in phase space y]; however, since a detailed explanation about the beats has not 

■ ■ ■ ' been presented until the present time, we judge that they deserve a deeper investigation. 

Recently, Chountasis and Vourdas |5|] showed that the Weyl function is an important tool for quantum interference 
effects. In particular, they studied the Wigner and Weyl functions for a superposition of to quantum states \si), where 
each function is decomposed into diagonal and nondiagonal terms, where the nondiagonal term is responsible for the 
interference between the states \si). Adopting a different approach and using the formalism developed in 0, and in 
order to shed more light on the origin of the beats predicted in Q , here we decompose the squeezed state HE into 
three functions: the two marginal (one for p and the other for q) and the correlation distribution functions (MDE 
and CDF, respectively). Our results corroborate the phase-space interference concept and complement the graphical 
treatment proposed by Mandal (t) . 
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This paper is organized as follows. In section II we discuss the solution of the differential pseudo-diffusion equation 
(see reference g), and in section III we define the phase space MDF and CDF. The formal and numerical results 
are given in section IV, where we show that oscillations and beats are present in the CDF function. Section V 
is dedicated to a summary and conclusions. Two appendices are also presented, containing calculational details. 
Appendix A contains the main steps to calculate the MDFs by direct integration, and in appendix B we calculate the 
CDF. 

II. THE PSEUDO-DIFFUSION EQUATION AND ITS SOLUTIONS 

The mapping of the statistical operator p (describing a state of the EM field) on the squeezed-states representation 
(SSR) permits us to write the HF P as follows [g: 

P{p, q; A, (b) - Tr [pn(p, q; ()] - Tr [p,U{Pr, Qr] A)] = F,(p,, g,; A) , (1) 

where 

n(p,(z;C) = |m;C)(to;CI -R(0/2)n(p„g,;y)Rt(0/2) (2) 

is a projection operator, R((/)/2) = exp I -^ a^a) is the rotation operator, qr = qcos{(j)/2) + psm{(j)/2) and pr =^ 

pcos{(j)/2) — qsm{(j)/2) are the rotated phase space variables expressed in terms of the old ones, p^ = IV((j)/2) p R(0/2) 
and A = e~^^ (0 < A < oo). Now, if one considers the mixed state p — X^^o^" !"')("• I' diagonal in the Fock basis, 
Pj. will be invariant under rotations since R(^/2)|7i) = e^'^^/^jn). Consequently, the associated HF is given by 
P{p,q;X,(j)) = P{pr,qr',X)- This relation is useful in the sense that it is sufficient to consider the calculation of the 
unrotated distribution P{p,q;X), with variables changed from {p,q) to {pr,qr) in the final result, respectively. The 
number state |?t.)('^| is a typical example where this relation can be used. In Q we demonstrated that P(p, g; A, (?i>) 
satisfies the partial differential equation 

T{p,q-\4>)P{p,q-\,^)^Q, (3) 

where 

r(p, g; A, ,/>) = —- ^ I [A^ cos2(0/2) - sin2(0/2)] — + [A^ sin2(0/2) - cos2(0/2)] — 

-(^' + 1)^-^5^} (4) 

is a linear differential operator. For (/) = 0, equation (S) is similar to the diffusion equation in two dimensions where the 
parameter A plays the role of time. In this situation, since the diffusion coefficients have opposite signs, the equation 
describes a diffusive (infusive) process in the p (g) variable. For this reason, it has been called the pseudo- diffusion 
equation H,^. 

Here we consider the formal solution of equation (g|), obtained by the Fourier transform (FT) method, written as 
an integral equation with the kernel depending on the squeeze and rotation parameters, 

P(p,g;A,0) = r ^ e^(''P-««)i^(C,77;A,0)P(e,r;) 



OO 



'^^^'^ J(...-«,.)^(^^^.^^0)P.(^,r;) 



, 2^ 

Pr{pr,qr-A) ■ (5) 



The kernel is 



K{^,7j;X,cb) = exp L^ {[Asin2(0/2) - cos^Wm^ + [Acos2(0/2) - sm\cf>/2)]rj^ + (A + l)sin0 ^rj}] (6) 

and P(f , 77) is the FT of the HF P{p, g) for A = 1 and = 0. Wc notice that K{^, rj; A, (p) is an unbounded function, 
responsible for the squeezing 'propagation' of an 'initial' function P(j), q) to P(p, g; A, (f) for any values of A and (j) in 



their domain. Thus the existence of a 'propagated' P{p, q; A, (p) depends on the functional form of P{£,, rj), since the 
integral in the first line in ^ for P{S,, il) = constant does not exist. From the pseudo-diffusion equation ^ and the 
linear differential operator F, equation (1q) shows the following symmetry properties 



and 



P{p, q; X, 4>) - P{q, -p; A, </) ± tt) = P{q,p; X-\ 



K{^, 77; A, 0) = K{f^, -e; A, </) ± tt) = K{r,, C; X-\ -^ 



(7) 



(8) 



Now, our aim is to show that the Glauber-Sudarshan distribution P'^{p, q; A, (f>) and HF P{p, q; A, (j)) are related by 



P^iP, q; A, 4') = A(p, q; A, 0)P(p, q; A, <?i) = P{p, q; -A, < 



with 

which can also be written as 
A(p, q; A, 0) = exp 



A(p, g; A, (j)) = exp 



Va^ + a-^ 

2 V dpy dq^ 



(9) 



(10) 



1 
2A 



since 



[A2 cos2(0/2) + sin2(0/2)]^ + [A^ sin2(0/2) + cos\cf>/2)]-^ 



5 = --'(^/2)^ + cos^(<^/2)|, - sin^^ , 



(A2-l)sin( 



dpdq 



g=cos^(0/2)^+sin^(0/2)^+sin< 



92 



dqdp 



Applying the differential operator A on (||), we get 

d£,d-q 



A(p, g; A, 0)P(p, g; A, </>) = 



27r 



A(p, g; A, ( 



Jivp^-iq^) 



2n 



K{^,7j;X,0)Pr{^,ij) 
K{^,rj;X,0)Pr{tv) 



K(i,V-X,0) 



d£,dr] 

P{p,q;-\ 



,^(^Vr-m.)K[^^^--X,Q)Pr{i,ri) 



(11) 



The second equality is obtained using the following relation 



A(p, g; A, < 



jiVPr-iQr) 



E 

fc=0 



2fefc! 



A"i— +A— ] e*(''P'-««') 
dq^ dpi 



00 _. 

2^ li 
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fc! 



^ (A-^e^ + X^f 



n fe 



ji.VPr-iqr) 



^ e5(^"'«'+^''') e*(''P'-"«9'-) 



and, from the definition of the kernel K{S^, rj; A, 0), we conclude that 



62 



(^^'i'+^r,')K{^,7r,X,0) = e3(^+i)('''+^"'«') = if(^,,7;-A,0) 



Thus, the DF P'^{p, g; A, 0) is directly obtained by changing the signal of the squeezing parameter A -^ —X in the HF. 
Consequently, this result shows that Pnip, q,X, (j)) does not exist as a bounded function for the number states, however 
it exists as an ultradistribution. Equation (||) is a generalization of previous results obtained in |^,0 for = 0. 



III. MARGINAL AND CORRELATION DISTRIBUTION FUNCTIONS 



The HF P{p, q; A, 0) can be written as a sum of two terms: the first is the product of the two MDF, in q and p, 
and describes the noncorrelated part; the second term is the CDF and contains the phase-space correlations Q. So, 
the HF @) can be written as 



where 



and 



Pip, q; A, 0) ^ R{p; A, (^)Q{q; A, (b) + C{p, q; A, 

dp 



Qiq;\ 



Rip; A, 



/27r 



dq 



P{p, q\ A, < 



Pip, q; A, . 



(12) 



(13) 



(14) 



are the MDFs, and C(p, q; X, 4>) is the CDF. However, the calculation of (|l^) and (|TJ) by direct integration displays 
difficulties when (j) ^Q (see appendix A). Thus, the aim of this section is to show that expressions for the MDFs can 
be obtained in a much simpler way by using the formalism of section II. 

Substituting the right hand side (RHS) of the first line of equation (||) into the integrands of equations (|l^) and 
(|l4|), and then carrying out the integrations we get 



and 



where 



and 



Qiq;X,(t>) 



Rip; A, < 



d^ 



drj 

^2^ 



e-'^'>kQi^;\,4>)QiO 



e'^nRip;\,<jy)Rip) , 



fcQ(^;A,0) =exp|-^ [Asin2(0/2)-cos2(</)/2)]r 



fcfl(?7;A, 0) = exp 



A-1 
4A 



[Acos2(0/2)-sin2(0/2)]?7^ 



(15) 



(16) 



(17) 



(18) 



are the reduced kernels responsible for the 'propagation' of the squeezing. The functions Qi^) and Rirf) are the 
respective FT of the Husimi functions Qiq) and Rip) for A = 1 (absence of squeezing). Futhermore, the equations 
([15) and (|l^) arc solutions of the partial differential equations 



d A2sin2(0/2)-cos2(0/2) d^ 
dX 4A2 g^ 

d A2cos2(0/2)-sin2(<^/2) d^ 



Q(g;A,0) = O 
i?(p;A,0) = O 



(19) 
(20) 



d\ 4A2 9p2 

In analogy to (|^) , the reduced kernels kq and kn have the symmetry properties 

kQ{R)ix;X,(j)) = fcfl(Q)(a;;A"\0) = fci^(Q) (x; A, tt ± 0) , (21) 

which reflect directly into the MDFs, 

Qiq; A, 0) = Riq; A-\ 0) = Riq; A, ^ ± </)) , 

Rip; A, (f>) = Qip; X-\ 4>) = Qip; A, ^ ± 0) . (22) 

Consequently, the calculation of Qiq; A, </>) is sufficient for determining the function Rip; A, (p), and vice- versa. 



Now we analyze the structure of the kernel (||) , which can be factorized as 



(23) 



where the first two factors on the RHS (reduced kernels) 'propagate' the initial HF in an independent way, i.e., if in 
the 'initial' (A = 1) HF the phase-space variables are not correlated, they will remain as such for any other value of 
A. The factor 



kc{Lm\(t>) ==cxp 



A^-1 



4A 



■ sin (/) 1 ^rj 



(24) 



introduces additional (or new) correlations into an 'initial' HF when <j) ^ nvr, with n G IN. Otherwise, we obtain 
kc = 1 and K[£^, t]; A, m:) — fcQ(^; A, mr)kfj{ri; A, m:), one reduced kernel for each variable. 

As a consequence of the factorization (p3) it is interesting to rewrite the CDF C{p, q; A, 0) as a sum of two terms. 



C{p, q- A, 0) = C(l) (p, q- A, 0) + C(2) (p, g; A, , 



defined as 



and 



C(i)(p,g;A,^) = 



2tt 



MnP~i'i)K{i,r^-\4>)C{tv) 



C<^'\p,q;X,^) = 



^ ekvP-ii)kQ(^- A, 0)fc^(^; A, 0) [kcii, 77; A, </)) - 1] R{v)Q{0 ■ 



(25) 



(26) 



(27) 



Here we assume that the HF contains 'initial' correlations [see equation (^] with its FT being P(^, 77) = R{'>])Q{S,) + 
C{^,r]), and that 'propagation' of correlations originates from two sources. The first, the RHS of equation (p6), is 
responsible for the 'propagation' of squeezing into the 'initial' correlations C{^,ri). In the second, equation (2^), 
'propagation' occurs only for 4> ^ mr when additional correlations are created into the 'initial' FT of the uncorrelated 
part of the HF, R{r])Q{S^). In appendix B, the correlations C'"'^' and C^^' are obtained for the Fock states in the SSR. 



IV. FOCK STATES IN THE SQUEEZED STATES REPRESENTATION 

The density operator p„ = l'^)('^l niapped in the coherent states representation yields a Poisson distribution |ll| 

1 /p2 _|_ q2 



Pn{p,q)^\{pq\n)\' 



il 



exp 



p2 _|_ q2 



with n = 0, 1, 2, . . .. So, the respective Fourier transform 

T7 ' 



^"(^'^) = W^E 



fc=0 



W^^-l^. ""'' (ti) ^''-'' \W '"^ ("^ 



(28) 



(29) 



where 7im(x) is the Hermite polynomial, represents an initial step for calculating the Husimi function Pn{p, q; X, (j)) 
in the squeezed states representation. In fact, substituting equation (E9|) into (M) and evaluating the integrations over 
C andr; O, §7.374-8], we get 



Pnip,q;X,, 



2%/A /A-l\" 1 



A + 1 V A + 1 / 2"?i! 



7Y„ 



Xqr + ipr 



exp 



Mr + pI 

A + 1 



(30) 



where qr and Pr are the rotated variables defined in section II. This expression was inittialy obtained in M], and 
later used by Schleich et al. ||l| in the oscillatory behavior study of the distribution P„(p, (7;A,0). Now, using the 
mathematical relation |l4| 



|H„(z)|2 = 2"n!^ {-lfLi-"^\2x') 4lf )(-2y2) (^ ^ x + iy) 



k=0 



in which Cm (x) is the associated Laguerre polynomial, equation (BG) can be written in an equivalent form, 



Pn{p,q;\4>) = 



2VA /A-1 



A+1 V A+1 



n n 



k=0 



A2-1 



C, 



(-1/2) 



2p^ 



A2 



exp 



Ml + pI 
A + 1 



The Husimi function Qnio) is obtained with the help of equation (|28|), i.e., 



whose Fourier transform is given by 



dp 



Pnip,q) =exp 



ry2k 



9 / Z^^n-k \^J) 2fefc| 



,(-1/2) 



k=0 



exp 



2 



fc=0 

Substituting this result into equation (|l5| ) and doing the integration with respect to ^, we get 

' " A - 1 cos2(0/2) - A sin2(0/2) 



Q„(q;A,0) = 



2A 



.777:;^E 



xC 



(A + l)[cos2(0/2) + Asin2(0/2)] ^^ 

(-1/2) \ 2A2g2 

fe 



(-i)'4-f^(o) 



(A2 - l)[cos4(0/2) - A2 sin*((/./2)] 



exp 



A + lcos2(0/2) + Asin2(0/2) 

V 

{X + l)[cos2(0/2) + Asin^(0/2)] 



In order to obtain the MDF Rnip', A, 
Rn{p;\(j)) = 



2A 



we only need the symmetry properties (P2), 

"A-lAcos2(0/2)-sin2(0/2) 



E^ 



x£. 



(A + l)[Acos2(</>/2)+sin2(0/2)] f^^ 
2A2p2 



(-1/2) 



(0) 



(-1/2) 



(A2 - 1)[A2 cos4(0/2) - sin^(0/2)] 



exp 



A + lAcos2(^/2) + sin2(</)/2) 

V 

(A + l)[Acos2(0/2) + sin2(0/2)] 



(31) 



(32) 



(33) 



(34) 



(35) 



Figures l(a)-(d) show the three-dimensional plots of Pn{p, q; X, (f>) versus n and (p, for A — 21,201,1/21,1/201, 
respectively. The plane = in figure 1 (a) corresponds to the oscillations pointed out in M , which depend strongly 
on (j), showing a periodicity of n. Now, for A = 201 [figure 1(b)] we observe the occurence of rich structures, although 
the beats pointed out in |4| can not be perceived. In fact, they are revealed in figures 2 and 3. Figures 2(a)-(f) show 
the plots of Pn versus n for — 85°, ...,90° and A = 201, where the beat structure becomes evident; however, it 
disappears at angles close to 90°. Following the arguments presented in 0] and corroborated by Mandal m, this beat 
structure is a consequence of the quantum interference in phase space. Figures 3(a)-(f) show the plots of C„ versus 
n for the same parameters used in figures 2(a)-(f), where the beat structure is present again. This fact connects the 
correlations and interference effects in phase space, and provides further insights to the phenomenon. Moreover, we 
observe a similar kind of plots for A = 1/21 and A = 1/201, figures 1(c) and 1(d), respectively, where now they are 
shifted by tt/2. 



V. SUMMARY AND CONCLUSIONS 

We have considered the Husimi function P{p, q; A, 4>) with emphasis on the marginal and correlation distribution 
functions, showing that all three satisfy the pseudo-diffusion equations if one considers that the squeezing parameter A 
plays the role of a time. The solution, obtained from the Fourier transform method, permits calculating P{p, q; A, 0), 
given an 'initial' Husimi function P{p, q) with a kernel Kjj, ry; A, 4>) responsible by the propagation of squeezing. The 
decomposition of the kernel in three factors, equation (|23|), permits writing the CDF as a sum of two terms, having 



different interpretations: the first term, equation (uq), is the propagation of 'initial' correlations contained in P(p, q); 
whereas the second term, equation (P7|), is responsible for introducing additional correlations into the uncorrelated 
'initial' product of the MDFs Q{q)R{p). 

Finally we remind that the formal procedure employed throughout this paper is advantageous if compared with the 
direct and lengthy calculation exposed in appendix A. In the specific case of the number state, the decomposition of 
the CDF into two terms should permit to investigate more thouroughly the origin of beats. Although having attained 
a formal expression for both (see appendix B), the numerical calculation presents difhculties due to its complexity. 

Multimode-squeezed-states representation can also be considered within the present formalism. In particular, M. 
Selvadoray et al. |lj] studied the two-mode-squeezed-state photon distribution. Again they verified the presence of 
beats. In this case, the correlation distribution function plays a crucial role in the understanding of this effect. 
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APPENDIX A: MDFS FOR THE NUMBER STATES BY DIRECT INTEGRATION 

The usual procedure to calculate the MDFs for the number states consists in the integration of equation 
respect the variables p or q, respectively, 



Q„(q;A, 0) 
Rn{p; X, (f>) 



dp 



2tt 
dq 



Pn{p,q;X,(t)) , 



P„(p,g;A,0) 



OD with 

(Al) 
(A2) 



We calculate ( |AlD by direct integration and present some properties inherent to the marginal distributions. 
Consider initially the integral representation of the Hcrmite polynomial |16[| 



7i„(z) = (-2i)"e/ / ^w"e 



du 



n „ — u +2luz 



which permits to write (l3(l) in a more convenient form: 



, ,, 2VA (-2)" /A-l\" 
F„(p, g; A, </>) = ■^— ;- , ( .;— ;- ) exp 



sin2(,/)/2)-Acos2(</>/2) 



A+1 n\ V^ + 1, 

cos2((/)/2)- Asin2((/)/2) 



X exp 



A-1 



A + 1 



P 



A-1 
('uw)"expJ -(u^ +v^) + 



A-1 



sm (p pq 



°° dudv , ,„ r . 2 2\ 2 



^/A^^T 



[{u - v) sm{(j)/2) + iX{u + v) cos(0/2)] q 



X exp 



VA^^T 



[{u — v) cos{(j)/2) — iX{u + v) sm{(f>/2)] p 



Substituting (A2) into ( |Al| ) and integrating with respect to p, we obtain 



g„(g;A,0) = V^i-p(^^)"e-("^)^ 



a/tt 



'•°° dx „ 
X / —= X exp ■ 



a; _ i ( /3g + i _ (^1 „ ^ ) y 



for values of the squeeze parameter in the intervals 1 < A < cot^((/)/2) or cot^(0/2) < A < 1, where 

and (i — 



X 



X 



Y (A + 1) [cos2(</./2) + Asin2(0/2)] 
Integration over the variable x leads to |1J, §3.462-4] 



Q„(g;A,0) = V2^i-|^(^^y'e-("^)^ 
Now, using relation |lj| 



(A - 1) [cos2(0/2) - Asin2(0/2) 









-{V-Iaqf 



(A3) 



(A4) 



(A5) 



Uniz^w) ^J2^n-kiO) {2wr-^nk{z) 



fc=0 



for the Hermite polynomial present in (A_5), 



Hn 



. P /. a- 



we get 



^^ + ^2^1^ /^: 



^2a2 /a 



Ec-.(o) 



A;=0 






■HkiPq) 



Q4.;A,0) = ^(^^j e-(-^^EC-^.(o) 



(3 fa 



2a \P^ 



71 — k 



Hk{l3q)'H2n-k{aq) 



(A6) 



Note that (|Aq ) can be written in a compact form and equivalent to equation (|34|), i.e., in terms of the associated 
Laguerre polynomial. For this purpose it is necessary to verify the equality 



22" 



-| '"■ /0\/-i \ iL — a iL 



s=0 



/?; \P' 



fc=0 



and to determine the coefficients Cnkia_J3). 

The relation established by Bailey |l^ for the product of Hermite polynomials, 



(A7) 






7 2-^1 ( TO + 1, — /; 771 — 2j + 1; 



a2 + 52 



X 'Hm+i-2] (^Va2 + 62 xj 



where 2-^1(01,02; 03; 2) is the hypergeometric function, permits to verify the equality (A7) through 

k\ 



'U(R\'U f \ (-l)"(2n-.)! ^ 

ns{f3q)n2n-s{aq) = 2n-s^^p-j2n-s 1^ 



fc=0 



(2fc-s)!(n-/s)! 



(2Q,)2fc^2(«-fc) 



1 



X 2.F1 ( 2n - ,s + 1, -,s; 2fc - ,s + 1; ^ ) 4 ^/'^ [2(a/3g)21 



and to determine the coefficients c„fc(a,/3), 

(-1)" fc! (2a)2'=/32("-fc) 



c„fe(a,/?) 



22" {n-ky. (2a2/?2)" ^ s!(n - s)!(2fc - s)! 



E 



(2n-s)! 



2«)^M^-i 



X 2.?^i ( 2n - s + 1, -s; 2/c - s + 1, 



1 

2/32 



(A8) 



(A9) 



In fact, the sum present in (A9) can be performed since we used the following relation for the hypergeometric function 



1 \ s!(2fc-s)! 
2J^i 2n-s + l,-s;2/e-s + l;— - ' - ^ ' 



E(-i) 



; {2n~s + ly. (2^2) 



'2\~l 



2/32; (2n-s)! ^' ' [s - iy.{2k ^ s + ly U 



(AlO) 



Then, substituting (AlO) into ([A9|) and calculating the sums, we obtain a simple expression for the coefficients 



c„.(a,/3) = (-l) [2(n~fc)]!! UJ = ^^'^ ^-^ (^\/3 



2fc / \ 2/c 

fcr(-i/2)/n^ f " 



(All) 



So, t he marginal distribution ( [A6| ) can be expressed such as equation (p3) . Adopting analogous procedure for equation 
(A2) we obtain (p5|), however the values of squeeze parameter are restricted in the intervals tan'^{(j)/2) > A > 1 or 
ta[?(0/2) < A < 1. 



In addition to the symmetry relations established in section III [see equation (22)], the MDFs exhibit the following 
properties: 



{i) lim Qn{q;X,0)^\^n{q)r and lim i?„(p; A,0) = |$„(p)|2 , 

A^oo A— >0 



E 

n=0 



(m) VQ„(q;A, 



^i?„(p;A,^) = V/:ti/^)(0) 



n=0 



n=0 



(A12) 
(A13) 



The first property does not characterize the Husimi function as a probabihty distribution but only emphasizes the 
character of quasiprobabihty distribution: hmit vahies of A recover the squared nioduh of wavefunctions |lO,n8[. 
With respect to the second property, it is a direct consequence of the scalar product for the squeezed states, i.e., 
^^^g P„(p, g; A, 0) = 1. In fact, the divergence is a consequence of the integration step of this relation over the 
variables p or q. Now, considering the normalization of equation (pOJ), /^^ -|^ Pnip, 9; A, 0) — 1, we obtain the third 
property for MPDFs: 



(Hi) 



dq 

v/2^ 



Qn{q\\4') = 1 



and 



dp 



i?„(p;A,0) = l. 



(A14) 



APPENDIX B: THE COMPONENTS OF THE CORRELATION DISTRIBUTION FUNCTION 

The components of the correlation distribution function C(p, g; A, </>) can also be written as 

C(l) {p, q- A, 4>) = P{P, q; A, 0) - C(3) (p, q; A, 0) , 

C(2) (p, g; A, ^) - C(3) (p, g; A, 0) - Q(g; A, 0)i?(p; A, ^) , 

with 



C(3)(p,g;A,^) = 



2n 



K{i,r^-X,4>)Q(^-l)R{in-l) 



(Bl) 

(B2) 



(B3) 



Present in both equations (Bl) and (B2), the term C'^^-' is responsible for introducing correlations into the marginal 
distributions product through the relation 



C(3)(p,g;A,</>)=exp 
where we used the result: 

For = 0, equations (B1)-(B3) simplify to 



A^-1 
4A 



92 



4A 



■ sin ( 



dpdq 



dpdq 



Q{q;X,^)R{p;X,(f>) 



e^(''P-«'')fcQ(e;A,0)fc«('7;A,0). 



(B4) 



C(i)(p,9;A,0)=C(p,g;A,0), 

C(2)(p,9;A,0)-0, 

C(3)(p,g;A,0)=g(q;A,0)i?(p;A,0). 

In this appendix we calculate the component C^^' for the number states and, consequently, the components C'^' and 
C'2' can be totally determined. 

In order to simplify the calculations, let us initially consider the marginal and correlation distribution functions 
expressed in terms of the parameters ai, a2 and 0:3: 



n 

Qn{q;X,<t>) = y2^e-("^^)^^ 4^:^,(0) (-2a?)'= 4-1/^) [{a.qf] , 

n 

R^{p; A, </.) = ^ e-(-^)^ ^ ci%{0) {-2alf ci'"^^ [{a,pf] , 

00 / p.21 
C^^\p,q;X,^)^J2^—^Q4q;X,^)R„(p;X,cf>), 



(B5) 
(B6) 

(B7) 



with 



ai = 



(A + l)[cos2(r/)/2)+Asin2(</)/2)] ' 



a2 



as = 



(A + l)[Acos2(,/,/2) + sin2(0/2)] ' 



1a" 



SllK 



and OiJ+ a\ = 2afa|. Expressions (B5) and ( jBg ) are the alternative way of writing the marginal distributions ( |34| ) 
and (|3|), respectively, since we used the properties of the associated Laguerre polynomials. Moreover, using the 
relation 

those expressions permit to calculate the component Ci , given the result: 



Cf (p,g;A,0) = 2aia, e-[(-«)^+("-)^l J] 4^1,(0) ^^ Y. ^""^^ 



fc=0 



a|^ 
2fefc! 



2™m! 



X > — I ^ Hi+2k{aiq)Tii+27n{a2p) 






(B8) 



The infinity sum present in (B8) can be done by using the formula |l9, §5.12.2.1] 



oo j-k 



t 1 



f x~2ty 



Atxy- 4*2 (a;2 + v"^) 



-, {m.n} 



Hn- 



l-4f2 
y — 2tx 



^ r!4^U(0) 4:2,(0) (4ir 



r=0 



where {m, n} stands for the minor of m and n, which leads to 

2Q!ia2 



|i|<l/2, 



C(3)(p,g;A,0) 



^1 - (2Q;ia2a3)^ 



exp 



afq ~ Aa'ia'^azpq + a2P 



2^2 



fc=0 



a? 



1 — (2aia2a3)^ 



2™to! 



1 — (2aia2a3) 



1 - (2Q;ia2a3)2 

m {2fc,2m} 

5] r!4t,(0)42-.(0)(4aia2a3) 



/ aig-2aia^a3p A / a2P - 2a^a2a3g 

X lt2k-7- — ;, lt2m-r 



a/1 - (2aia2a3)^ 



^/l - {2aia2a3y 



(B9) 
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FIG. 1. Plots of P„(p,q;X,(f>) vs. n and (p with p^ + q^ = 98. Figures l(a)-(d) correspond to A = 21,201, 1/21 and 1/201, 
respectively. 

FIG. 2. Plots of Pn{p,q; \,<j}) vs. n with <j} — 85°,..., 90° and A = 201. The phase-space variables are transformed into 
q = 7v2cos^ and p — 7\/2sin6', which leads to qr = 7\/2cos(^ — 4'/'^) ^nd pr — 7-\/2sin(S — 4>/2). For mapping the parameters 
used by Dutta et al., we fixed 6 — 3(j>/2. 

FIG. 3. Plots of Cn(p, q; A, <j}) vs. n for the same parameters set used in figure 2, where we see the presence of beats again. 
This fact corroborates the phase-space interference concept since correlations and interference effects are closely connected. 
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